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Abstract 

The charmless nonleptonic B — ► r\K decay is an useful probe to test the strong interaction 
dynamics part of the Standard Model. Within the QCD factorization framework, we analyse 
this particular decay by using the most recent data of BR(B + — > i]K + ) and A cp (B + — ► r)K + ) 
available in the Heavy Flavour Averaging Group (HFAG) webcite. Using these data we 
constraint the unitarity angle 7 and 77 — 77' mixing angle 0. We find that such constraint is 
scale dependent, e.g. for 7 = 70°, the data for BR(B+ -> 77^+) (= (2.6±0.5) x 1(T 6 ) suggests 
that 9 should lie in between —46° and —44°, —30° and —26° for /_t = mb/2. For /1 = nib, 
the same BR data (with the same 7) suggests that should lie in between —54° and —50°, 
—30° and —26°. The allowed region followed from A cp — — 0.25± 0.14, is found to be rather 
wider than that obtained from the BR(B+ -> r)K + ) data. For 7 = 70° and 6 = -21.3°, we 
find BR(B~(+) -> r)K-(+)) = 2.93(4.91) x 1CT 6 and A cp = -0.252 at fi = m b , and about 
2.98(5.72) x 10~ 6 and —0.315 corresponding to /1 = mb/2. We investigate the role of the 
power corrections in such constraints. 
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1 Intrduction 



The nonleptonic weak hadronic decays of the B meson is an useful probe to test the Standard Model(SM), 
particularly the dynamics of it's strong interaction (QCD) part. Due to the non-perturbative features, 
the amplitudes for such hadronic S-decays, are difficult to calculate directly from the QCD lagrangian. 
For the phenomenological study, several factorization hypotheses, like naive- factorization(NF), general- 
ized factorization(GF) and QCD factorization(QCDF) were introduced and has been quite successful in 
explaining the data for exclusive B decays, like B — ► PP, PV, VV, PT, VT (where P, V and T stands 
for pseudoscalar, vector and tensor meson) reported by the CLEO, BABAR and BELLE collaborations 
0121 El- Out °f all these hypothesis, the QCD factorization hypothesis 0] particularly, has become very 
much popular in recent years because of it's successful explanation of the several B factories data. Study- 
ing the nonleptonic _B~(+) — > rjK~( + ) decay within the QCD factorization, exploring it's CP-violation 
aspects, is a well motivated topic and is the main concern of the present work. 

The only source of the CP violation in the Standard Model(SM) is the CKM matrix Vckm 0- which 
is arised due to the misalignment of the mass and weak interaction eigenstates. It is expressed via the 
charge current-current interaction lagrangian Cf^, defined by 

C frS = —J= ( U L C L t L ) J^VcKM 

In the Wolfenstein parametrization Vckm is given by 

1-±A 2 




v ud 


v us 


Vub 


Vcd 


V cs 


v cb 


v td 


Vts 


Vti, 




-X AX 2 I +C(A 4 ), (2) 

AX 3 (1 - p - iif) 

where A = V us — sin8 c , 9 C , the Cabibbo mixing angle. The angle 7 related to the phase of the CKM 
element V u b i.e. 7 = arg(V* b ) (which follows from V ub = \V u b\e~ 11 — AX 3 (p — if])), amounts to CP 
violation in the SM. The angle 7 is one of the three angles a, (3 and 7 of the unitarity traingle and 
any nonzero value of it (resulting into non zero area of the unitarity traingle) yields the CP violation in 
the SM. Hence, a precise measurement of the unitarity angle 7 can be possible from the CP violating 
transition, which is being considered in the present work. There exist bounds on the angle 7 and one 
such bound is followed from the time-dependent CP asymmetry measurement. The current average 
experimental value of the SV-n- which is about = —0.49 ± 0.27, yields [Hj 

7 =(66111)" or (1743)° , (3) 

where both set of solutions are consistent with each other and the first set is consistent with the unitarity- 
traingle fit. Our next concern is the mixing angle 9 between the psedoscalar 77 and 77' mesons, a matter 
of great interst from the time when the SU (3) flavour symmetry was proposed. In the simplest scenario 
where the 77 and 77' mesons do not mix with other pseudoscalar mesons, like excited quarkonium states, 
gluonium or exotics, the 77, 77' wave functions can be written as 

1 77) = cosO \r) S ) - sinO \r] ), (4) 
1 77') = sinO 1 77s) + cos6 1 770), (5) 

where the SU(3) basis states |tjs) and 1 770} are given by 

1 778) = —=\uu + dd — 2ss), (6) 
V6 

1 770) = —^=\uu + dd + ss). (7) 

From the accumulated data of the decay widths r[?7 -> 77] = (0.46±0.04) KeV, T[r)' -> 77] = (4.26±0.19) 
KeV and r[7r° -► 77] = (7.7 ± 0.55) eV, one finds |7|, 

9= -21.3° ±2.5°. (8) 
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In the present analysis we will treat the unitarity angle 7 and mixing angle 9 as the unknown parameters 
and use the most recent data for the direct CP asymmetry A cp (B + — ► i]K + ) and BR(B + — > 7]K + ) [H] 
to put constraints on them. 

The organization of the paper is as follows. Section [21 begins with the general discussion of the 
nonleptonic decay of a B meson within the QCD factorization framework, which includes the discussion 
of several non-factorizable corrections, e.g. vertex, penguin, hard-spectator and the weak annihilation 
corrections. They are present in the QCD factorization framework, but absent in the simplest naive 
factorization(NF) framework. We obtain the SM decay amplitude Msm{B + — > r)K + ), BR(B + — ► r]K + ) 
and A cp (B + — > rjK + ). Section |31 is fully devoted to the numerical analysis. After describing several 
numerical inputs i.e. the CKM matrix elements, effective coefficients, quarks masss, decay constants, 
form factors, we discuss the constraints in the 9 — 7 contour plane, which is being obtained by using 
the most recent data of BR(B + — ► t]K + ) and A cp (B + — ► i]K + ), available in the HFAG webcite [HJ- In 
Section^ we summarize our results and made our conclusion. 



2 B 7) K decay within the QCDF framework: 

The most general effective weak Hamiltonian Tt^f =1 for the non-leptonic AB — 1 transitions can be 
expressed via the operator product expansion (OPE) ^HI 
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v U bV* s (doy + c 2 o%) + v cb v* (ciOf + C2 o 2 c ) - v tb v t * 8 ]T 



H.c. (9) 



Here c^s are the wilson coefficients and the 4-quarks current-current, gluonic and electroweak penguin 
operators are defined by 

• current-current operators: 

O" = {ub)v~A(su) V -A O2 = {u a bp)v-A{spu a )v-A, 

0\ = (cb) V -A(sc)v~A O2 = {c a bf})v-A{sf}C a )v-A; (10) 

• QCD-penguin operators: 

O3 = (sb)v-A^2(qq)v-A, 4 = (s a bp) V ^A^2(q/3<la)v-A, 

q q 

O5 = {sb) v _ A ^2(qq)v+A, 6 = {'s a b [j ) v ^ A ^{q l3 q a ) v+A , (11) 



q q 



• electroweak-penguin operators: 

O7 = -(sb) V -A^2e q (qq) v+A , 8 = -(s a bp) V -A e q (q p q a ) v+A , 
q q 

°9 = -^{~sb)v-A^e q (ciq) V -A, O w = -{s a bp)v~A^2e q ^f}q a )v~A, (12) 
q q 

where a, (3 are the SU(3) color indices, V ± A correspond to 7 M (1 ± 7 5 ) and the wilson coefficients c,'s, 
in which the QCD correction to weak interaction is encoded, are evaluated at the scale 0([i ~ mb). e 
and g are respectively QED and QCD coupling constants and T a 's are 5/7(3) color matrices. For the 
penguin operators, O3, . . . , O10, the sum over q runs over different quark flavors, active at /1 ~ mb, i.e. 
g e{u, d, s, c, b}. Note that the gluonic penguin operators C*3_6 contributes largely to B —> r/K. The 
contributions coming from the electroweak penguin operators are not so significant and will be neglected 
in the present analysis. 

The wilson coefficients CjS, quite well-known at the NLO order accuracy, is evaluated at [i ~ mt> and 
is within the perturbative control. The nontriviality (amounting uncertainity) exists in the evaluation of 
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the hadronic matrix elements of the operators in the effective Hamiltonian comprising B — > r\K transi- 
tion. Several approximations are made in order to have a control over this. Within naive factorization 
framework, one assumes the absence of the order 0(a s ) QCD correction to the hadronic matrix element 
and the validity of working in the heavy quark limit (i.e. rrib 3> Liqcd)- Within the naive factorization, 
the hadronic matrix of a B meson decay to a pair of mesons Mi and Mi can be factorized and be written 
as the product of two quarks bilinear currents j\ and ji 

{MxM 2 \Oi\~B) NF = {Mx\jx\B) x (MaljalO) (13) 

which gives rise the transition form factor and the decay constant, respectively. The operators O^, 
constituting the effective Hamiltonian, is specific to a particular B decay. In the present case they are 
represented by Eqns. i|10lllll2|) . Corrections to the naive factorization arises from the "non-factorizable" 
0(a s ) QCD correction to the hadronic matrix element by means of vertex, penguin, hard-spectator 
and as well as the weak annihilation contributions. The "QCD factorization" framework, which can 
be looked as the order 0(a s ) corrected version of the naive factorization framework, encodes all these 
0(a s ) corrections to the hadronic matrix element (|13fl . Within QCD factorization, the above two body 
nonleptonic decay amplitude (Eqn. I|13|) . in the heavy quark limit, can be generalized as 0] 

{MiMi\Oi\B)QaDP = y"Ff^ Ml (m|) / duT^u) $ M2 (u) + (M x <-► M 2 ) 

Jo 



1 



+ / dCdudvT l u ^,u,v)^ B (O^Mdv)^M 2 {u) (14) 
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where Mi and Mi corresponds to the light mesons. The factor F B ^ M i( M 2) corresponds to the B — > Mi 
transition form factor, while the <E>x corresponds to the Light Cone Distribution Amplitudes (LCDA) 
for the meson X for the quark-antiquark Fock states. The perturbatively calculable quantities T 1 ^ 11 in 
Eqn. Q14f> . which are known as the hard-scattering kernels, contain essentially the 0(a s ) QCD correction 
to the matrix element. which starts at 0(a®), corresponds to the vertex, penguin correction, while 
T^ 1 which starts at 0{a\) corresponds to the hard-spectator correction which arises due to the exchange 
of a hard gluon between the spectator quark inside the B meson with the quark of the emitted Mi meson. 
Note that the weak annihilation contribution are not shown here, we will present them later seperately. 
Now no QCD correction implies the vanishing of the second term of Eqn. i|14|) . while turns into a 
constant. After performing the relevant integration (Eqn. (jl^ l with the leading twist-2 LCDA's, one in 
the heavy quark limit ends up a expression, which much looks like: form factor x decay constnt, agrees 
with the naive factorization result (Eqn. lfl3|) 1. Following the above discussion, the generic B — > MiM 2 
decay amplitude, in the heavy quark limt mt, ^> Aqcd within the QCDF framwork, can be written as 



(M 1 M 2 \O i \B) QC DF = (MiMilO^B) 



NF 



m b 



(15) 



Note that in the heavy quark limt (m& 3> Aqcd) and at order 0(a s ), although the naive factorization 
framework breaks down (due to the presence of the second term), one can still calculates the the cor- 
rections systemetically by finding the corresponding corrections in the short-distance coefficients and the 
LCDA's of the mesons. 

In the QCDF framework, the amplitude for the B — * r\K can be expressed as 

M S m(B -> V K) = M f SM (B -> rjK) + M a SM (B -> V K), (16) 

where 

M f SM (B^ v K) = ^ X> aP i {vK\Oi\B) N F, (17) 

v p—u,c i—1 
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i=6 



M a SM (s^ V K) = -^} B f Kfri Y J Y. l P h 



(18) 



In above l p — V p bV* s with V's being the CKM matrix elements, Gf, the Fermi decay constant, while 
f r (r = B, K,rj,rj') stands for the meson decay constants. The non-factorizable vertex, penguin and 
hard spectator corrections are encoded in the effective coefficients aF i which appeared in the first term of 
the amplitude ( Eqn. JIEJ). The effective coefficients aP i at the next-to-leading order in a s {\i) takes the 
following form (for a generic nonleptonic B — ► M\M 2 decay) [H] 



af(Mi,M 2 ) = ct(ji) + 



Ci±l( 



(A*) \ . c i± i(/x) C F a s (n) 



-in 



At: 2 

V i (M 2 ) + —H l (M 1 ,M 2 ) 



Pf(M 2 ), (19) 



where Ci{fj), the wilson coefficients with i(= 1, ...6) is odd(even), the upper(lower) sign is applied and 
for the current-current operator the superscript p is to be dropped (to avoid confusion) . Here N c = 3 is 
the colour factor, Gf = the quadratic Casimir Invariant for SU(N) (N — N c = 3 here) and Cj's 

are the wilson coefficients. The quantities Vi, Hi and Pf respectively stands for the "non-factorizable" 
vertex, hard spectator and penguin corrections and their complete expressions can be found in 
The hard spectator functions reads as 



Hi{Kn) 

for i = 1-4, 



is /if/," 



$b(0 dx dy 



xy 



xy 



(20) 



(0)/," 



dx / dy 



*v(x)* K (y) %(x)<t>f(y) 



for i = 5 with and 77 i (M 1 Af 2 ) = for i = 6. In above j = 1 



xy 



xy 



(21) 



, y) and the term 



2ml 



m i ,(/i)(m,+m 



-y^y, (where m q corresponds to the average of the up and down quark masses) is the chiral 



enhancement factor. The Q>x($x ) corresponds to twist-2 (twist-3) LCDA's of the meson X. The weak 
annihilation contribution's b^s corresponding to Eqn. (|18l) . can be expressed as (following @],[n|) 

Gf ,i , p Gf 



b 2 



N ^A h bl = 



c 3 A\ + c 5 {Al + A{) + N c c 6 A f 3 



(22) 
(23) 



where 61, b 2 corresponds to the current-current annihilation, while 63, b\, the gluon penguin annihilation. 
Here A^ is the factorizable annihilation amplitude which is being induced from the (S — P) x (S + P) 
type operators, whereas the non-factorizable annihilation amplitudes A\ 2 3 are induced from the (V — 
A) X (V — A), (V + A) x (V + A) and (S— P) x (S + P) type operators. Their explicit expressions, which 
can be found in 0],|ni, reads as 



4 

4 

A' 



4 



Al = 



4 



0. 

0, 



dxdy i §M 2 { x )®MAy) 



dxdy |$ A/2 (x)$ Ml (2/) 



dxdy 



1 



2/(1 - xy) 



— 9 



x(l — xy) xy 2 



2[J 



+ rf rf $2(x)$^(y)£- 



2 



] 



:ry 



,(3) , 



2x 



7ra.(/i) / dxdy ^ rf 1 $ Ma (1) (</) 



,(3) 



2(1 + 5) 



x 2 y 



+ 



rf * Ml (y)*$Ax) 



2(1 + 2/) 



xy 



(24) 
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for a generic B — > M\Mi decay. In our case of interest, the corresponding expressions are obtained by 
setting Mi = K(rj) and Mi = rj(K). Note that the power suppressed twist-3 LCDA $^ of the hard 
spectator and annihilation contributions suffers from the end point divergence Xh,a = Jq — > which can 
be phenomenolgically parametrized as [S] 



X 



H 



V \ l> h 



(25) 



where Aft ~ 0.5 GeV, pu,A are free parameters expected to be about of order ph,a — 0(1) and 
4>h,a e [0, 2tt}. With all these QCD corrections in hand, we are now ready to give the 0(a s ) corrected 
expressions for the B — ► r]K decay amplitude (i.e. Eqn. ifTTjl) and within QCDF they reads (for B = B~) 
as jZj 



M f SM (B- 



rjK-) = 



^§ \Vu b v: s 



°2 + a l~ 



F 



«) f. 



K 



m 



K 



^0 



-v tb v* 



2a p 3 - 2og 



K) ft 



V 2 



V cb V c * 8 a^ 



fC 

' fu 
J r) 



P 2 



m s (mb - m s ) 



fs 

f% m s (mb-m s ) 



+ 



2a\m 1 K 



(m b 



Tor — mz 



K 



F, 



(m 2 K ) f. 



F 



(K-\sb_\B-) (r)\uu-\0) 



(26) 



where i>_ (and u_ ) corresponds to 7^(1 — 7s)& (and 7 M (1 — 7s)u). The transition form factor and decay 
constants are defined as 



(K-(p')\sb-\B-(p)) = {K-{p')\s^b\B-{p)) 



(p+pT- 



nr. 



m 



K 



q»F^ K tf) 



(27) 



(Tj( g )|««_|0) = -(r y (g)|S7 AI 7 5U |0) = i/^. (28) 

where ms, ttik corresponds to B,K meons masses and f™, the r\ meson decay constant. F^i* K (l 2 ) 
are the B — > K transition form factors evaluated at q 2 = ml. Note that in the effective coefficients 
a p (Eqn. (|19fl ) appearing in the above amplitude (Eqn. 1261) ) contains the non-factorizable vertex (V*) 
and penguin (Pf) corrections, which are evaluated at fi = mb (rrib/2) and the hard spectator (Hi) and 
annihiation corrections (given by Eqn. 1)22(1 ). evaluated at fi = ph ~ VT^h> where L/j ~ 0.5 GeV [Hj- It 
is now straightforward to write down the branching ratio (BR) for the B~ — > r\K~ decay as 



BR(B- - ryX-) = \M(B~ - 7,A-)| a il>!)) 



where A4(B — ► ^AT ) is the decay amplitude, tb, the B meson life time and p c , the cm. momentum 
of the r\ and K~ mesons in the B~ rest frame, is given by 



_ V ( m B - ( m v ~ m K -) 2 ) (m 2 B - (m, + m K -) 2 ) 
Pc — -z ■ [ov) 

ZniB 

The BR's BR(B + — ► r/K + ) is obtained from BR(B~ — > r)K~) by performing the relevant changes under 
CP transformations. The CP-asymmetry A cp , according to the standard convention, is defined as 



T(B~ -> rjK-) - T{B+ -> ^+) 

r(#- tja:-) + r(s+ t]A:+) ' 



4* - ; T;„ + _^ - (3i) 
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Table 1: Data for the B + — ► r/K + decay mode. The error bars are at la limit; the upper limit 
at 90% CL.For A cv we follow the standard convention as mentioned in the text. 



Final State 


BR x 10 e 


-A-cp 


r]K + 


2.6 ±0.5 


-0.25 ±0.14 



3 Numerical Analysis 

The BR and CP-assymmetry A cp data for the B + — > ij K + decay are available in the HFAG05 website 
|H| and is given in Tabled The decay amplitudes depend on the effective coefficients a,'s, CKM matrix 
elements, quark masses and the non-perturbative inputs like hadronic form factors, decay constants etc. 

3.1 CKM matrix elements 

The parameters A,l,p and rj of the CKM matrix Vckm i n the Wolfenstein parametrization are set 
at the following values. We employ A and I = sm6 c at the values of 0.815 and 0.2205 in our analysis. 
The other parameters are found to be p = \/~p 2 + rj 2 cosj and r\ — \J~p 2 + rj 2 sin-f with 7 = cirg(V* b ) (an 
unknown parameter in our analysis) and y/'p 2 + rj 2 = 0.3854 9s . Here p = p(l — 4p) and rj = rj(l — 4j-) 

EH- 

3.2 Effective coefficients a^, quark masses, decay constants and form factors 

The NLO wilson coefficients in the NDR scheme, obtained in the paper by Beneke et al. [H],are cataloged in 
TableEl The effective coefficients a^s (i = 1,2,..., 10) in Ean. (|26|l . can be obtained by using Ean. l|19j) from 

Table 2: AB = 1 wilson coefficients at /i = ^{m b ) ~ 2.1(4.2) GeV for m t = 170 GeV, a = 1/129 
and A^L = 225 MeV in the NDR scheme EQ. 



NLO 


c\ 


C2 


C3 


c 4 


C5 


C6 


pi = m b /2 


1.137 


-0.295 


0.021 


-0.051 


0.010 


-0.065 


pL = m b 


1.081 


-0.190 


0.014 


-0.036 


0.009 


-0.042 



these NLO wilson coefficients. The wilson coeffficients Ci(p) and coupling constant a s {p) (in Eqn. I|19|l) 
are evaluated at p — m^ira^jl). The scale chosen for the vertex, penguin, hard-spectator and annihilation 
term are described above (see the discussion before Eqn. (|29|l b 

For the quark's constituent masses, we use m& = 4.2 GeV, m c = 1.5 GeV and m s — 0.50 GeV and 
m u = rrid ~ 0.2 GeV, the scale independent quantities [HllZj, which appear in the loop integral and for the 
current masses, we use their scale dependent values as listed in Table which appear in the factorized 

Table 3: The scale dependent current quark masses which are taken from 0)0- 





m b (p) 


m s (p) 


m u {p) 


m d {p) 


p = m b /2 


4.88 


0.122 


0.0413 x m,(ji) 


0.0413 x m s (p) 


pL = m b 


4.20 


0.090 


0.0413 x m s {p) 


0.0413 x m s (p) 
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amplitude after making the use of equation of motion of quarks. 
The decay constants (in GeV) are given by (|H],|Zj) 

f B = 0.20, f n = 0.131, f K = 0.160, f°, = 0.0058, f c q = 0.00093 (32) 

and we use fo = 1.10 f n and fg = 1.34/^ with the one angle mixing scheme (in n — n' sector) to obtain 

n /s cosQ f sin9 



= h cos9 _ f sine 

if] m rx 1 J n 



V6 



V3 



V6 



V3 



(33) 



with e, the mixing angle, another unknown parameter in our analysis. 

We use the following non-perturbative form factors values in our analysis [H], 



(0) = 0.28, F V (0) = 0.34 



F, 



(0) = F — (0) 



cosO sine 



For the distribution amplitudes (LCDA), we use the asymptotic form for the pseudoscalar mesons 



(34) 



$(3) 
V 



6x(l -x), $k = 6x(l - x), twist - 2 LCDA 



(3) 
K 



1, 



twist - 3 LCDA. 



The B meson wave function 



$ B {p) = N B p 2 {l- pfexp 



1 / pm B 



2 V ujb 

with the normalization constant Nb is being determined from the condition 



f $ fl (p)dp = l, 
Jo 



(35) 
(36) 



(37) 



(38) 



with m B = 5.278 GeV and uj b = 0.25 GeV Q31- 
3.3 Results and Discussions 

With all the numerical inputs in hand, we are now ready to discuss our results. Working within the QCD 
factorization framework, we made use of the BR(B + — > nK + ) and CP-asymmetry A cp (B + — > nK + ) |S] 
data to obtain constraints in the — 7 plane. As we will see that the bound depends crucially on the scale 
[i at which different nonfactorizable 0{a s ) corrections are evaluated. The impact of the nonfactorizable 
weak annihilation correction, which is usually power suppressed, will also be considered while obtaining 
such constraints. 



3.3.1 Constraints in the 6 — 7 plane: with (without) the annihilation terms 

We have seen that the effective coefficients en's fEan ll9(l . contains the wilson coefficients a's, coupling 
constant a s {p) and several nonfactorizable corrections, which crucially depends on the renormalization 
scale For the wilson coefficients Ci in <Zj, we use their NLO values for different p which are displayed 
in Table [3 The current quark mass which depends on the scale p., are displayed in Table |3J For the 
running coupling a s (/j,), we use the oneloop expression for a s (p) as [TU] 



4?r 



l3 Log p 2 /A. 



MS 



(39) 



where /3o = (^N c — |n/) with N c — 3, the colour and the typical QCD scale parameter Aj^(Aqcd) — 
226 MeV with n/ = 5, the number of active quark flavour at /1 = nib and similarly K-jj^[Kqcd) = 372 
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MeV with rif = 4 for u. — nib 1 2 JHj- Using these as inputs, we now use the BR and A cp data ( Table 
^1 to obtain constraint in the 8 and 7 plane. Before to obtain such constraints, let us see how the CP 
asymmetry A cp varies with 7 or 8. In Figures l(a,b) and 2(c,d), we made such plots. From Figures l(a,b) 
( standing respectively for /1 — m&/2 and nib with the mixing angle 8 = —21.3° [7|), we see that within 
the ±1(7 deviation with respect to the the central value Ajif" = —0.25, 7 is allowed to varry in between 
+20° and +100° (one range), 140° and +170° (other range). The negative 7 is ruled out at the ±ler 
level as far as the CP asymmetry data is concerned. From Figures 2(c,d), in which the unitarity angle 

7 is chosen as +70° the scale /z, is being set at m&/2 and nib, respectively, we find that at the ±ler 
crror away from the central value A^ 1 , the mixing angle 8 is allowed to lie in between —45° and —40° 
(one range), —25° and —10° (other range). However, the A cp data allows another region, in which 8 is 
positive. However, the constraints on 8, which follows from the r\ — ► 77, 7/ — > 77 and 7r° — > 77 decay 
studies, suggest that 8 should be negaive i.e. 

8= (-21.3 ±2.5)°. (40) 

In order to have the overall consistency with other's analysis, we will not explore the positive 8 region. 
From the above analysis, where one angle is varied, while the other fixed, we know their sizes and signs. 
We will now explore the general possibility in which both 8 and 7 is allowed to vary. In this case, by 
using the A cp (B + — > r/K + ) data, as well as the BR(B + — ► r/K + ) data, it is possible to constrain the two 
angles 8 and 7 simultaneously by finding the contour plots in the 8 — 7 plane. 

In Figures 3(e, f), we obtain such contour plots corresponding to fi = nib/2 and m&. The central 
curve corresponds to the region allowed by the A cp — —0.25 (i.e. the central value), while the next(outer) 
to the central curve corresponds to A cp = — 0.11(= —0.25 + 0.14) and the region between the outer and 
central curves, is allowed at the +la level. The innermost curve corresponds to A cp = —0.39 and the 
region between this curve and the central one, is allowed at the — la level. Note that a different set of u, 
choice results into the change in the shape and hence the area of the allowed regions as a whole. At the 
±1(7 level, in either case, quite a large range of both 8 and 7's are allowed: e.g. while the mixing angle 

8 varies in between —42° to —6°, the unitarity angle 7 can changes largely in between, say, e.g. 10° to 
175° at the ±1<7 level. As for example, say 7 = 70°, the contour plot suggest that 8 should lie in between 
—44° and —40° (left zone), —24° and —6° (right zone) for the scale /x = rrib/2, while in between— 42° and 
— 40°(left zone), —26° and —8° (right zone) for the scale /i = m?,. 

In Figures 4(g, h), the 8 - 7 contour plots corresponds to BR(B+ -> nK+) = (2.6 ±0.5) x 10~ 6 . The 
renormalization scale fi is chosen at rrib/2 (Figure 4g) and nib (Figure 4h). There are two regions in this 
plot. For each region, the central curve corresponds to the central value BR(B + — » nK + ) = 2.6 x 10~ 6 , 
while the next (outer i.e. rightwards for the right side, and leftwards in the left side) to the central curve 
corresponds to BR(B + — > r]K + ) = 3.1 x 10~ 6 and the region between the outer and central curves (on 
either side), is allowed at the +la level. The two innermost curves corresponds to BR(B + — > nK + ) = 
2.1 x 10 -6 and the region between this curve and the central one (in either side), is allowed at the — la 
level. The allowed contour region followed from the BR data, is rather narrow than that obtained from 
the A cp data. Within ±lcr error bar, while 7 can change widely in between 0° to 180°, 8 changes very 
little: it change only by ±1° with respect to the central value (corresponding to the central curve). Note 
that a different /1 shift the 8 — 7 allowed contour region (see Figure 4g and 4h). As an example, say 
7 = 70°, the contour plot suggest that 8 should lie in between —46° and —44° (left region) and in between 
—30° and —26° (right region) for the scale fi = nib/2, while for fi = nib, it lies in between —54° and —50° 
(left region) and —30° and —26° (right region). 

To see the explicit scale dependence of the BR and CP asymmetry Acp results, let us calculate them 
for some chosen u. values. At u, — nib/2, with 8 = —21.3° and 7 = 70°, a point well within the allowed 
region, we find BR{B~^ -> r?iC"W) = 2.98(5.72) x 10" 6 and A cp (B- -> r}K~) = -0.315 when the 
annihilation term is taken into consideration. If the annihilation term is switched off, the correponding 
predictions turns out to be BR{B~^ -> nK'^) = 3.0(5.83) x 10~ 6 (a change about 2%(11%)) and 
A cp {B~ — > r)K~) = —0.319 (a change of ~ 4%). Similarly, for fj, = nib, with the same 7 and 8 value 
and considering the the annihilation contribution, we find BR(B-M -► nK-t+y) = 2.93(4.91) x 10~ 6 
and A cp (B- -> nR-) = -0.252. Without the annihilation term, we find BR(B~ (+ ^ -> jyJf-W) = 
2.94(4.93) x 10~ 6 (a change about 1%(2%)) and A cp {B~ -> nK~) = -0.253 (a change of - 1%). So, the 
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impact of power correction (annihilation contribution) on BR and A cp at different fj, is not so significant 
and consequently, in constraining the parameter space of 9 and 7. 

In Figures 5(i, j) and 6(k, 1), we have plotted A cp and BR{B + — » rjK + ) as a function of 9 and 7 with 
different set of the renormalization scale i.e. /i = mi,/ 2 (left one) and m& (right one). Note the differences 
arises due to the choice in renormalization scale: in Figures 5(i,j) larger negative A cp may arise at /i = mi, 
than that due to \x = m&/2 and in Figures 6(k,l), in the same way large BR also arise at /1 = mi,/ 2 in 
comparison to fi — mi,. Finally, from the Figures 5(i, j), we see that a large negative A cp , consistent with 
the present data, allows 9 to lie around —40°, for a wide range of 7, consistent with our finding. 

4 Summary and Conclusion 

We have investgated the charmless nonleptonic B — > r/K decay within the QCD factorization framework 
in the Standard Model. Two crucial things in the decay amplitudes are: (i) the unitarity angle 7, and 
(ii) the Tj — rj' mixing angle 9, which possesses several independent constraints followed from different 
studies. By knowing either 9 or 7 from others analysis and using the HFAG A cp data, one can get some 
idea about the allowed range of 7 or 9. We considered the most general scenario in which both 9 and 
7 are being treated as the unknown parameters. We use the BR(B + — > i]K + ) and A cp (B + — > rjK + ) 
(available in the HFAG webcite) to constrain the 9 — 7 contour plane. Using the A cp data, at the ±lcr 
level, we found that when the angle 9 varies in between —42° to —6° (in two ranges), the angle 7 can vary 
in between 0° to 180°. For the BR(B+ -> r)K+) = (2.6 ±0.5) x 1CT 6 , with 7 = 70°, we found that 9 can 
vary in between —46° and —44°, —30° and —26° for /i — mb/2. For /j, = mi,, the corresponding allowed 
range is in between —54° and —50°, —30° and —26° with the same choice of 7. We find the diference 
in the allowed contour region of 9 which follows from the A cp and BR constraints, a crucial finding of 
the present study. We made our analysis for different set of renormalization scales i.e. /i's, which results 
into change in the allowed region in the 9 — 7 contour plane. This is an important observation of this 
work. The effect of power correction correction is also investigated and found to be not so significant 
in constraining the 9 — 7 contour plane. This is so because it's (power correction) impact on A cp and 
BR(B-M _> riR-M) arc about 1% and 1%(2%), respectively. 
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gamma 



(a) (b) 

Figure l[a, b]. Plot showing the CP asymmetry A cp (B + — » r/K + ) as a function of the unitarity angle 
gamma (7, given in degree). The most recent data J^j for A cp — —0.25 ± 0.14. The Figure la and lb 
corresponds to the following choice of the renormalization scale: /1 = ^ (Fig. la) and [i = (Fig. lb). 
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(c) (d) 

Figure 2[c, d]. Plot showing the CP asymmetry A cp (B + — > r/K + ) as a function of the rj-rj' mixing 
angle theta (9, given in degree). The CP-asymmetry A cp = —0.25 ±0.14 accoriding to the HFAG webcite 
fSy. The Figure 2c and 2d corresponds to the choice of the renormalization scale: /i — ^ (Fig. 2c) and 
/! = rrif, (Fig. 2d). 
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(e) (f) 

Figure 3[e, f]. Contour plots in the theta (0) - gamma (•y) plane (angles are given in degree) using 
A cp = —0.25 ± 0.14 with the renormalization scale /i = rrib/2 (Fig.Se) and m& (Fig.Sf), respectively. 
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Figure 4[g, h]. Contour plots in the theta (9) - gamma fa) plane (angles are given in degree) us- 
ing BR(B + — > rjK + ) — (2.6 ± 0.5) x 10 -6 with the renormalization scale fi = mb/2 (Fig.^g) and 
rrib(Fig.J h h), respectively. 
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Figure 5[i, j]. Showing the CP-asymmetry A cp (B + — > rjK + ) as a function of the mixing angle theta 
(9) and the unitarity angle gammafj). The angles are given in degree. The left one corresponds to the 
renormalization scale /i = mb/2 (Fig.5i), whereas the right one, fi = n%b (Fig.Sj). 
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Figure 6[k, 1]. Showing the BR(B + — > r)K + ) as a function of the mixing angle theta (degree) and the 
unitarity angle gamma (degree). The angles are in degree. The left one corresponds to the renormalization 
scale fi — rrib/2 (Fig. 6k), whereas, the right one, /i = m& (Fig. 61). 
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